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ABSTRACT 

Context. In the time-ordered data (TOD) files of the Wilkinson Microwave Anisotropy Probe (WMAP) observations of 
the cosmic microwave background (CMB), there is an undocumented timing offset of —25.6 ms between the spacecraft 
attitude and radio flux density timestamps in the Meta and Science Data Tables, respectively. If the offset induced an 
^ I error during calibration of the raw TOD, then estimates of the WMAP CMB quadrupole might be substantially in 

error. 

Aims. A timing error during calibration would not only induce an artificial quadrupole-like signal in the mean sky map, 
it would also add variance per pixel. This variance would be present in the calibrated TOD. Low-resolution map- making 

1/^ ' as a function of timing offset should show a minimum variance for the correct timing offset. 

Methods. Three years of the calibrated, filtered WMAP 3- year TOD are compiled into sky maps at HEALPix resolution 

A^side ~ 8, individually for each of the K, Ka, Q, V and W band differencing assemblies (DA's), as a function of timing 

offset. The median per map of the temperature fluctuation variance per pixel is calculated and minimised against timing 

offset. 

Results. Minima are clearly present. The timing offsets that minimise the median variance are —38 ± 8 ms (K, Ka), 

^ . -27±3 ms (Q), -43±8 ms (V), and -47±194 ms (W), i.e. an average of -30±3 ms, where the WMAP collaboration's 

Mh' preferred offset is ± 1.7 ms. A non-parametric bootstrap analysis rejects the latter at a significance of 99.999%. The 

hypothesis of a —25.6 ms offset, suggested by Liu, Xiong & Li from the TOD file timing offset, is consistent with these 
minima. 

Conclusions. It is difficult to avoid the conclusion that the WMAP calibrated TOD and inferred maps are wrongly 
calibrated. CMB quadrupole estimates {3/tt)C2 based on the incorrectly calibrated TOD are overestimated by roughly 
64 ±6% (KQ85 mask) to 94 ± 10% (KQ75 mask). Ideally, the WMAP map- making pipelines should be redone starting 
from the uncalibrated TOD and using the —25.6 ms timing offset correction. 

> 

O 

en ' 1. Introduction the WMAP Explanatory Supplement (Section 3.1, Limon 

'/"J ' „^,„,^^ 6t al. I2010p does not refer to the —25.6 ms offset. 

r I . Liu & Li (|2010ap made sky maps out of the time-ordered 

(^ ■ data (TOD) files of the Wilkinson Microwave Anisotropy Liu et al. (|2010aD presented a simple model for the ef- 

Q I Probe (WMAP) observations of t he cos mic microwave feet of a timing error on subtraction of the Doppler-induced 

'background (CMB) (Bennett et al. I2003bp . following the dipole signal, which is celestial-position-dependent because 

procedure recommended by the WMAP collaboration, it is a dipole and time-dependent because of the space- 

They found a considerably weaker CMB quadrupole sig- craft's orbit around the Sun, demonstrating that an artifi- 

^ . nal than that estimated by the WMAP collaboration, but cial qu adrupole-like signal would be generated. Moss et al. 

;h ' were not aware of the reason for the difference, given that (|2010p presented a toy model that confirmed that a timing 

5^ ] their software pipehne satisfied many tests. Later, Liu et al. offset of -25.6 ms would induce an artificial quadrupole ap- 

(|2010ap found that (i) they had used a different choice of proximately aligned with what had been previously consid- 

timing interpolation than that of the WMAP collabora- ered to be the cosmological CMB quadru pole sig nal. There 

tion, and (u) there is a timing offset of -25.6 ms between was n o clear consensus between Liu et al. (j2010ap and Moss 

the spacecraft attitude and radio fiux density timestamps, et al. (j2010p regarding the fraction of the quadrupole that 

recorded in the Meta and Science Data Tables, respectively, would be artificial, 
in the TOD filesQ Individual observations in the Q, V and 

W bands last for 102.4 ms, 76.8 ms, and 51.2 ms, respec- ^^ the timing offset induced an error in compiling the cal- 

tively (Section 3.2, Limon et al. [20T0ll . so the offset corre- i^rated TOD mto maps, then this should mtroduce a blur- 

sponds to half of a W band observing interval or a quarter "^g effect at a few-aremmute scale. This effect would not 

of a Q band observing interval. The 21 April 2010 version of be easy to see by eye m sky maps. However, statistically, m 

full-sky maps created without ignoring bright objects, the 

hypothesis that a timing offset introduced an error in the 

^ The sign is chosen as the Science Data Table time minus the compilation of the calibrated TOD into m aps was excluded 

Meta Data Table time. to very high significance (Roukema I2010D . 
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Nevertheless, this method did not exclude the possibil- 
ity that a timing offset could have induced a quadrupole-like 
artefact during the calibration of the uncalibrated TOD. 
Indeed, Jarosik et al. (Section 2.4.1, Figs 3, 4, Jarosik et al. 
I2007P showed that an error in the gain model could lead to 
a quadrupole-like artefact in inferred sky maps. Although 
the Doppler dipole is not subtracted from the data during 
the calibration step, it is the signal used for the calibration 
itself. 

In principle, it would be possible to recalibrate the un- 
calibrated TOD using several different values of the timing 
offset, remake sky maps from the newly calibrated TOD, 
and test a statistic of these sky maps that shows which 
timing offset is optimal. However, a simpler approach is 
possible. 

A timing error during calibration would not only induce 
an artificial quadrupole-like signal in the mean. It would 
also add variance to the sky map signal per pixel. This 
would be present in the calibrated TOD. In other words, 
the effect of the error would be to add a small noise compo- 
nent to the signal per pixel. Since the calibrated TOD con- 
sist of individual observations, as do the uncalibrated TOD, 
reversing an incorrect timing offset when map-making from 
the calibrated TOD should approximately correct the dom- 
inant signal, the dipole, reducing this variance. Here, mini- 
malisation of the variance per pixel, as a function of timing 
offset, is carried out. Maps of both 5T and [ST)^ are gener- 
ated from the calibrated TOD as a function of timing offset 
in order that variance maps can be inferred. The timing off- 
set that gives the minimum variance should be that which 
is correct. 

The would-be timing-offset-induced calibration error 
and the way that an approximate reversal of the er- 
ror should reduce the per-pixel variance are explained in 
Sect. 12. ll A toy model illustrating the variance reduction is 
presented in Sect. 12.21 The observational analysis and vari- 
ance minimisation method are described in Sect. 12.31 The 
timing offset hypotheses are summarised quantitatively in 
Table [TJ Results are presented in Sect. [3J with the main 
results in Table [H A rough estimate of the relevance of 
the results for the CMB quadrupole is made in Sect. 13.11 
The discussion section deals with sensitivity to the fit- 
ting method (Sect. 14. 1[) . consistency between wavebands 
(Sect. 14. 2p . the possible relevance of sidelobe and other 
beam effects (Sect. 14. 3p . consistency with the results of us- 
ing an alternative method (Sect. 14. 4p . and various caveats 
(Sect. l43)) . Conclusions are given in Sect.[5l Gaussian error 
distributions are assumed throughout, except where other- 
wise stated. 



2. Method 

2.1. Model of would-be calibration error 

The full calibration procedure is described in Hinshaw et al. 
(|2003p . along with details for the 3-year data in Jarosik 
et al. (|2007p . Calibration to the dipole signal plays a funda- 
mental role in this procedure. The direction and amplitude 
of this signal vary slightly throughout the year, mainly as 
a result of the orbit of the WMAP spacecraft around the 
Solar System barycentre. 

The dipole signal used in the calibration procedure in- 
cludes both an estimate of the contribution from the Solar 
System barycentre's motion with respect to the comoving 



coordinate system, and a component due to the WMAP 
spacecraft's motion around the Solar System barycentre 
(Sect. 2.2.1 Hinshaw et al. 120031) . The latter is measured 
to high precision, presumably as a function of a time stan- 
dard, e.g. Julian days, consistent with that used in the main 
mapmaking procedure. 

Let us suppose that a timing offset is present when cal- 
culating the dipole that is used for the calibration of one 
hour of TOD, i.e. the times associated with the TOD space- 
craft attitude quaternions are used rather than the times 
associated with the TOD. Suppose that the only signal ob- 
served is that of the true dipole, and that the random errors 
in measurement are negligible. Over the approximately an- 
nular region corresponding to one hour's observations, the 
result of the calibration will be to shift the input signal in 
sky position to match the dipole calculated for a slightly 
incorrect position. That is, the gain and baseline will be 
wrongly estimated by a small amount that depends on the 
wrongly calculated dipole expected to be present in the ob- 
served annulus, in a way that on average corresponds to a 
positional shift. 

Over a full year, the set of observations in a given pixel 
(for some given pixel size) depends on the various hourly 
subsets of TOD of which that pixel is a member. The set 
of temperature fluctuation estimates in that pixel will have 
non-zero variance induced by the varying errors from differ- 
ent hourly calibrations. Thus, by calibrating with a timing 
offset, noise per pixel is added to observations that by as- 
sumption consist of the true dipole and negligible random 
measurement error. 

Now suppose that we have a set of wrongly calibrated 
observations (TOD), retaining the assumptions that the 
only observed signal is the dipole and that measurement 
error is negligible. Consider the meaning of calibration. In 
any hourly subset of TOD, the measured signal is offset and 
scaled in such a way that it matches the calibration signal. 
Thus, we know what the result of wrongly calibrating the 
input signal is. It is the assumed dipole. Hence, the wrongly 
calibrated values can be corrected towards the true val- 
ues by using knowledge of the scanning path, the assumed 
dipole, and the timing offset. This reduces the variance per 
pixel down towards the original, negligible measurement er- 
ror. The error cannot be reduced to zero, in this simplified 
case, because the calibration is not carried out on individ- 
ual observations (which would be meaningless), it is carried 
out on (usually) hourly subsets of the TOD. 

This can be expressed mathematically as follows. Given 
errorless observations 5T[(f){t),9{t)] of the dipole, the hourly 
calibration of these correct dipole observations to an incor- 
rectly timed dipole is equivalent to applying the functional 
gh ^ g, where g^ is the functional for the hourly calibration, 

g : dT[(j){t),e{t)] ^ ST[^{t - dt), e{t - St)] (1) 

is the functional for what would hypothetically consist of 
calibrating every observation individually, and the timing 
offset is St. This approximation should hold because the 
dipole is a smooth function with features only on large 
scales. Thus, the wrongly calibrated TOD can be written 
as ghiST). 

The inversion of the shift g^^ is 

g-^ : ST[(j){t),eit)] ^ ST[dp{t + St),e{t + St)]. 

Thus, 

g-^{g^[ST{(p,e)]] « g-^{g[ST{<j,M] 



(2) 
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Fig. 1. Minimisation of a median-per-map standard deviation 
per pixel a{5t) [Eq. (|16|l ] in a realisation of a toy simulation of 
would-be wrongly calibrated data with an input timing offset of 
St = 0.5, fit by a parabola over a range of timing offsets centred 
at St = 0.5. 



<5r(0,6i). 



(3) 



Hence, 

a{g-'{gi,[ST{cP,e)]}) 



« a[<5T((/.,0)] 


f/, ft 


< e 


< <j{g[ST{<l>,e 


)]} 



(4) 



where the variance a^ is defined for the set of observa- 
tions in a pixel at a fixed position {(f>,9). The relation 



a[5Tiq,,0)] 



< e is the assumption that the observa- 



tions have negligible (< e ^ 1 in any units) measurement 
error. Hence, given a set of wrongly calibrated observations 
6T' := gii{ST) of originally errorless observations 5T of the 
pure dipolc signal, minimisation of 



<j[9-\sr)] 



(5) 



as a function of St should yield the timing offset implicit in 
the wrong calibration. 

In reality, the signal is not expected to be a pure 
dipole signal, and the measurements are not errorless. 
Moreover, an iterative procedure involving both calibra- 
tion and mapmaking is used to optimise the calibration 
procedure (Sects 2.2, 2.3, Hinshaw ct al. I2003|) . and a more 
complex gain model is used in the calibration of the 3-year 
data (Sect 2.4, Jarosik et al. I2007[) . Thus, some additional 
error will be introduced into the above minimisation by the 
complex calibration procedure. 

Another potential complication is that use of g~^ will 
also induce a blurring effect at small (arcminute) scales 
fRoukema l20101 and references therein). At this small scale, 
application of g~^ should increase the variance per pixel if 
the dipole is assumed to be constant on this scale. Hence, in 
order to detect the original calibration error, a large enough 
pixel size must be used so that the dominant contribution 
to the variance per pixel is from positional offsets to the 
dipole rather than from blurring of point sources. 



2.2. Toy model 

A toy model illustrating the detection of a timing offset by 
minimising the variance per pixel can be set up as follows. 
Consider a fixed map vector m of A'pjx pixels, composed of 
a fixed sinusoidal "dipole" map (a one-dimensional simpli- 
fication of a dipole) defined 



"^o(j) '■— 3 sin 27r 



J 



N, 



pix 



(6) 



at pixel J and a component that includes CMB signal and 
noise, modelled (for simplicity) as Gaussian noise of zero 
mean and standard deviation 0.5, i.e. m is pixelwise defined 



"i(j):=mo(j)+0.5G(0,l). 



(7) 



The units can be thought of as mK. 

A TOD vector d representing A^tod successive ob- 
servations during one year can now be defined via an 
TVtod X TVn 



vpix mapping matrix M{5t) at timing offset 5t 



M{5t)m+ M{St) -MiiSty 



where • is the Hadamard produciQ, 



[Mi((5i)]„. := 0.3 sin 



27r- 



An 



■27r 



i + 5t 



Ntod 



(8) 



(9) 



is an A^pix x A^tod matrix representing the dipole compo- 
nent that varies during the year, and e is a column vector 
of size A^pix whose entries are 1, used for summation. The 
timing offset St may either be the correct value, in which 
case M{6t) is the correct mapping matrix and Mi{St) is the 
correct time- varying dipole, or St may be a guessed value. 
Let the mapping matrix M(St) be initialised to zero. 
For the i-th TOD observation, define the first beam column 
number 



JiiSt) 
mod 



70.5^^ 
360 P''' 



sin IGOtt 



i + St 
Ntod 



+ 



Nr 



TOD 



A^pix 



(10) 



where [.J is the floor function, set the second beam column 
number 



j2[jliSt)] 



mod 



.nm + 



141 
360^ 



rA^n 



,N„ 



and set matrix elements for the two beams 



M, 



iJliSt) 



^^^Mjiist)] 



l/(l + a;im) 
-1/(1 -a-i„,), 



(11) 



(12) 



where Xim = ±0.007 is a differencing assmbly imbalance 
parameter with a randomly chosen sign in a given sim- 
ulation. The motivation for these functions is to have a 
scanning pattern whose fractional sky coverage and frac- 
tional year coverage roughly correspond to those of the 
real observations, e.g. one beam varies approximately si- 
nusoidally in addition to a linear (S^) yearly orbit, and the 
second beam is separated from the first by a pixel distance 
of (141°/360°)Arpix. 



Entry- wise product. 
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The first and second moments of temperature fluctua- 
tion estimates calculated directly from the calibrated TOD 
d are maps (<5T}| , ^ and {{^T)^'^\ each of whose j-th 

pixel is 



Table 1. Hypotheses of timing offset 5t expressed in units of 
an observing interval At, and in ms. 



J2AM'^{St)U~M{5t)mn 



M{5t)-Mi{5tf e\] 



E.1 



(13) 



and 



E 



{St)i.d-M{St)mQ- M{St)- Mi{Stf ej] 



'g^,\ST)\<p,e 



iiSTn.giSt) 



{ST)\,JSt) ,(16) 



For each St, define the median per map of the standard 
deviation per pixel cr -i,xrr\ i n 



a{5t) :=^i/2(a^-^i(5j,)|^ J 



(17) 



Fitting a simple symmetrical function, i.e. a parabola, to 
<T{5t), should yield an estimate of 5t that minimises a{dt). 
This can be compared to the known value input to the 
toy model. Figure [T] shows the dependence of cr((5t) on 6t 
for a realisation of the model defined here, with A^pix — 
20, A^TOD = 2000, and an input timing offset 6t = 0.5. The 
range of timing offsets is chosen to be symmetrical around 
6t = 0.5. There is clearly a minimum in a close to St — 
0.5. Over 30 simulations with these same parameters, the 
median and standard error in the median of the St values 
that minimise the best-fit parabola to the median per map 
of <7{St) are 0.500 ± 0.003. 

2.3. Observational data and variance minimisation 



As in Roukema (pgTU| . only the three-year WMAP TOD 
are analysed, in order to reduce the computing load. 
However, the full-year, filtered, calibrated TOEJfl are anal- 
ysed (not just 198 or 199 days), for all three years. Analyses 
are made using the K, Ka, Q, V and W bands. Since there 
is only one differencing assembly (DA) for each of the K 
and Ka wavebands, giving only three one-year data sets 



^ |http : //lambda . gsf c . nasa . gov/product /map/dr2/ 1 
|tod_f cal_get . cfm] 



units : 
hypothesis 


band 


St an'" 
At At 


(St 


-0.5)Ai 
ms 


ms 


WMAP'' 


all 


0.5 -•= 




0.0 


1.7 



E.1 

(14) 
respectively, where the sums are taken over the rows i where 

{Al^{St) id - M{St) Too - M{St) ■ Mx{Stf ejV 7^ 0, 

(15) 

and (0, ff) is the celestial position of pixel j- Given the TOD 
d for an unknown timing offset, and the dependence of the 
mapping matrix on timing offset M(5t) from Eqs ([TU)) . ((TTj) . 
and (|T2|) . the standard deviation map can be estimated 
element-wise 



LXL2010 (i)'' K, Ka 


0.01 


-64.0 


1.7 


Q 


0.02 


-51.2 


1.7 


V 


0.02 


-38.4 


1.7 


w 


0.03 


-25.6 


1.7 


LXL2010 (ii)'= K, Ka 


0.30 0.01 


-25.6 


1.7 


Q 


0.25 0.02 


-25.6 


1.7 


V 


0.17 0.02 


-25.6 


1.7 


w 


0.00 0.03 


-25.6 


1.7 



^Bennett et al. (|2003a|l 's (Sect. 6.1.2) estimate of the "relative 

accuracy [that] can be achieved between the star tracker(s), gyro 

and the instrument" is adopted here for each of the hypotheses. 

''WMAP collaboration's preferred timing offset. 

'^Same per waveband as for Liu et al. H2010a|) hypotheses. 

''Liu et al. (|2010a|l 's hypothesis (i), i.e. described in terms of the 

beginning of an observing time interval. 

•^Liu et al. (|2010ap 's hypothesis (u), i.e. the TOD file timing 

offset. 



each, these are combined together to form a single set of 
six DA/year setsQ The Q and V bands each have two DAs, 
giving six DA/year sets each, and the W band has four 
DAs, giving twelve DA/year sets. 

For any DA/year set, for a range of timing offsets, the 
filtered, calibrated TOD are compiled using a patch to Liu 
et al. (|2010ap 's publicly available data analysis pipelined 
The patch is optimised for the present analysis and enables 
execution of the pipeline using the GNU Data Language 
(GDL)0 In contrast to the analysis in Roukema (j2010p . 
standard masking is retained in order to avoid planet ob- 
servations (daf_mask) and the Galactic PlaneQ The "pes- 
simistic" mode of downgrading the processing mask resolu- 
tion is used. In order to calculate variances, each mapmak- 
ing step creates both a mean signal map and a mean square 
signal map, as in Eqs (|13p and p4p . This is equivalent to 
using the first map iteration. 

Use of a highly-iterated map estimate would be equiv- 
alent to adding an extra term that would vary with St in 
a complex way, to these equations. This would add an ad- 
ditional source of noise to the statistic, decreasing its sta- 
tistical power. Nevertheless, given that the calibration is 
made using a complex iterative process both of sky maps 
and baseline and gain parameters, it should be useful to see 
if this gives a result that is statistically compatible with the 



■* For brevity, the K and Ka bands will sometimes be referred 
to below as a single band. 

^ |http : //cosmocof f ee . inf o/viewtopic . php?p=4525 [ 

http : //dpc . aire . org . cn/data/wmap/090 72731/release_vl/| 
source_code/vl/ and associated libraries 

http : //cosmo . torun . pl/GPLdownload/ 
^LLma pmaking_GDLpatches/LLmapmakiiig_GDLpat ches_0 . . 4 . tbz| 

' [h ttp : //lambda . gsf c . nasa . gov/dat a/map/ dr2/ 1 
|ancillary/wmap_processing_r9_mask_3yr_v2 . f its| 
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more accurate results using Eqs ([T^ and (fH)) . i.e. without 
addition of the extra term. Thus, an extra set of calcula- 
tions has been made using the 80-th mapmaking iteration 
for the bands with the smaller numbers of DA's and obser- 
vations. 

Liu et al. (|2010ap 's convention for labelling the timing 
offset, including the sign, is retained. That is, the center 
parameter for determining the timing offset is generalised 
to an arbitrary floating point value, written here as 6t, ex- 
pressed as a fraction of an observing time interval in a given 
band. Liu et al. (|2010a[) 's hypothesis of what is the correct 
offset, interpreted to mean that the true error occurred 
during the calibration step and not during the mapmak- 
ing step, includes two versions, labelled (i) and (ii) above 
(Sect.[T]). As (i) a "timing interpolation" , the hypothesis can 
be defined to mean St = in all bands. The WMAP collabo- 
ration's preferred timing offset corresponds to 5t = 0.5. The 
second version of Liu et al. (|2010ap 's hypothesis, (ii) above 
(Sect.d)), is the TOD file fixed timing offset relative to the 
WMAP collaboration's choice, i.e. {St-0.5)At = -25.6 ms, 
where At is the exposure time of a single observation in the 
K and Ka, Q, V, or W band, i.e. 128 ms, 102.4 ms, 76.8 ms, 
or 51.2 ms, respectively. Thus, in the W band, —25.6 ms 
relative to St ~ 0.5 corresponds to St ~ 0, while at lower 
frequencies, the offset fixed in milliseconds corresponds to 
values of St between and 0.5. The different hypotheses of 
the timing offset are summarised in Table [T] as a function 
of waveband. 

As explained in Sect. 12.11 the pixel sizes must be large 
enough to avoid introducing variance due to the incorrect 
inversion of small-scale signals (e.g. objects at the few- 
arcminute scale are blurred by g~^). Hence, a HEALPix 
(Gorski et al. ll999p resolution of A^sido = 8, i.e. pixel sizes of 
about 7.3°, is adopted here. This should reduce the chance 
that small-scale effects contaminate the effect of interest 
here. 

For each DA/year combination i in each band, maps 
of the signal {ST)\,g and the square signal ((<5T)^)| 

are calculated over —4 < St < 5, i.e. at values symmet- 
ric around St = 0.5, with smaller intervals closer to that 
value. This wide range in St is used because, as for the 
blurring effect analysed in Roukema (j2010[) . exaggerating 
the absolute timing offset should strengthen the effect of 
using a wrong value. 

For a given St, the two maps arc used to infer a map of 



the standard deviation per pixel a 



W(«T)U,e 



where ST now 



represents the calibrated (whether correctly or not) TOD. 
Let ai{St) be the median of this quantity over the map, 
as in Eq. ([T7)) . where invalid pixels (mostly in the Galactic 
Plane) are ignored. 

These "median standard deviations"|3 <Ti{St) are nor- 
malised over St internally within each sample i in order to 
give approximately equal weight to the different DA/year 
samples i in a given band. That is, the minimum Ai and 
maximum Bi of <Ti (St) are used to calculate 



^m := 



(7,{St) - Aj 
B, - A, 



(18) 



In order to estimate the minimum of Si(St) in a given 
DA/year combination in a given band, a smooth, sym- 



The standard deviation is calculated at a fixed position 
(0,6). The median of this quantity is calculated over all valid 
positions {4>,6). 
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Fig. 2. Mean and uncertainty of Si{5t) [cf. Eq. (fT8)l ] in the K 
and Ka, Q, V, and W bands, from top to bottom, respectively, as 
a function of 5t in units of observing intervals. For each band, 
symbols show the mean and standard error in the mean for 
the 6 or 12 DA/year combinations in that band. In addition 
to Si based on Oi defined as the median standard deviation in 
a map (0, labelled "p = 50%"), equivalent statistics based on 
the 40% (-I-) and 60% (x) percentiles are also shown. A least- 
squares best-fit parabola to the 50% curve is shown in each 
panel for convenience, but not used in estimating the minima 
Ti. The values 5t = Q and 5t — 0.5, favoured by Liu et al. 
(|2010a|) [version (i)] and the WMAP collaboration, respectively, 
are shown as vertical lines. Offsets in St of ±0.05 are applied to 
the 40% and 60% percentile statistics to reduce symbol overlap. 
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Table 2. Quadratic fit estimates of tfie timing offset [median 
IJ'i/2(ji) and standard error in tlie median o-jij/jC'^i)' ^1' GSl 
tliat minimises tlie standard deviation per pixel over 5t inter- 
vals symmetric around 5t = 0.5 in the range —4 < 5i < 5, 
for 6 DA/year combinations i in the K and Ka bands grouped 
together, for 6 DA/year combinations i in each of the Q and 
V bands separately, and for 12 DA/year combinations i in the 
W band, shown in units of an observing interval Ai, and in 
ms, compared to the WMAP collaboration preferred value of 
5t = 0.5. 



^1/2(1-0 0-Mi/2(^>) 


At ^1/2 (r, - 


-0.5)At a^,/,(..)At 


units : At At 


ms 


ms ms 


band 







C, Ka 


0.20 


0.06 128.0 


Q 


0.24 


0.03 104.2 


V 


-0.06 


0.11 76.8 


W 


-0.41 


3.74 52.1 



38.4 


8.3 


27.1 


2.8 


43.0 


8.1 


47.4 


194.8 



all 



0.22 



0.02 



-29.7 



2.5 



''In W, the best-fit parabola for one of the year/DA combina- 
tions had a maximum and was ignored. 

''Mean and standard error, weighted by the standard error, 
calculated either in observing interval units (ignoring the 
dependence of At on waveband), or in ms compared to (5t = 0.5. 



metric function is needed. The simplest obvious choice is 
a parabola. Since Si is normalised to the range [0,1], the 
parabola ai{5t)'^ -\-hi5t + Ci that least-squares best-fits Si{5t) 
is found, giving an estimate 



-h/{2ai 



(19) 



of the 6t value that minimises Si{5t). For a given band, each 
of the 6 (K and Ka, Q, or V) or 12 (W) DA/year combina- 
tions i yields an estimate r^. Under the assumption of statis- 
tical independence among the DA/year combinations, the 
median ^i/2{Ti) and standard error in the median u^j^ .^(Tijj 
give an estimate of the optimal timing offset for the inver- 
sion of a would-be timing offset during calibration of the 
uncalibrated TOD in that waveband. 

For the purposes of testing the sensitivity to the fit- 
ting method (Sect. 14. ip . additional calculations at 5t = 
—5, —2.5 were made in order to have a sample that is 
symmetric around 5t = 0, i.e. to test sampling sensitiv- 
ity, for all wavebands. Extra maps were also calculated for 
5t ~ —10, . . . , —6, 6, . . . , 10 in W to see if a wider range in 
5t gives a better defined minimum (Sect. 14. 2p . 

3. Results 

Calculations on 4-core, 2.4 GHz, 64-bit processors 
with 4 Gib RAM, using GDL-0.9~rc4 running under 
GNU/Linux, took about 2-6 hours per map, depending 
on wavebandEj Figure [2] clearly shows that the per pixel 



^ Gaussian error distributions are assumed, i.e. the standard 
error in the median is estimated as 1.253 times the standard 
error in the mean. 
'° The maps of the main 
of the a posteriori extra 



calculation and 
W calculations 
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Fig. 3. Histogram of the timing offset estimates fJ,i/2{Ti — 
0.5) At from the individual DA/year combinations, of which per- 
band statistics are listed in Table [2l and a Gaussian distribution 
appropriately normalised. 
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Fig. 4. Median standard deviation in /iK with zeropoints re- 
moved, i.e. ai{5t) — Ai as in Eq. H18|l . averaged (mean) over 
DA/year combinations in each of the K and Ka, Q, V, and W 
bands (symbols as indicated in the figure), as a function of tim- 
ing offset in ms relative to the WMAP collaboration hypothesis, 
i.e. ((5f — 0.5)At. Error bars are standard errors in the mean. The 
timing offset range is centred on that of the hypothesis favoured 
from Table m i.e. {St - 0.5) At = -25.6. 



variance of the maps is minimised near the preferred time 
offsets of the WMAP collaboration and Liu et al. (|2010ap . 
Time offsets of several observing intervals above or below 
these values clearly increase the per pixel variance. It is 
not clear in the figure whether any of the claimed offset 
hypotheses is significantly preferred as a sharp minimum. 
However, the approximate symmetry axes of the function 
shapes clearly lie at lower St than the central value 5t = 0.5. 
Tabic [2] shows the statistics of Tj, the minimum of Si for 
the individual DA/year combination for each waveband. 
Comparison with Table [T] shows that the WMAP collab- 
oration's preferred timing offset is rejected to high signifi- 
cance in the lower frequency bands, i.e. 4.5cr (K, Ka), 8.2cr 
(Q), and 5.2o- (V), and not constrained in the W band. 
Moreover, the K and Ka bands, and the Q band, also re- 
ject version (i) of Liu et al. (j2010ap 's hypothesis to high 
significance, i.e. 3. Oct (K, Ka) and 7.6cr (Q). Hence, both 



are available at http : //cosmo . t omn . pl/GPLdownload/ 1 
JLLmapmaking_GDLpatches/LXLof f set2_skymaps . tbz[ 
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the WMAP collaboration's preferred timing offset and ver- 
sion (i) of Liu et al. (j2010a|) 's hypothesis are rejected to 
high significance by minimisation of per pixel variance of 
maps made from the calibrated, filtered TOD. 

In contrast, the final two numerical columns in Tabled 
show that version (ii) of Liu et al. (|2010ap 's timing offset, 
defined as the TOD file timing offset of —25.6 ms, is con- 
sistent with all wavebands, at 1.5ct, 0.5cr, 2.1(T, and O.lcr in 
K and Ka, Q, V, and W, respectively. The weighted mean 
/^i/2(tj - 0.5)At = -29.7 ± 2.5 ms rejects the WMAP col- 
laboration's hypothesis at ll.Qcr and agrees with the fixed 
TOD file timing offset hypothesis at 1.6a. 

The above significance estimates assume that the statis- 
tic T is normally distributed. Figure [3] shows the 25 timing 
offset estimates, from individual 30 DA/year combinations, 
that lie in the range (—100, 100) ms. The distribution is not 
well-modelled by a Gaussian distribution. A Kolmogorov- 
Smirnov test with the mean and standard deviation of the 
data rejects the hypothesis that the distribution is Gaussian 
with p = 1.4%. Could this affect the significance in rejecting 
the St = 0.5 hypothesis? 

It is clear that it would be difficult to reconcile the dis- 
tribution in Fig. [3] with having a median of ms, especially 
given that the distribution is sampled many times. This 
can be tested more formally. To non-paramctrically test 
the hypothesis that the median of the distribution is ms, 
10^ bootstrap resamples from the distribution were made. 
These give a 99.999% confidence interval for the median of 
the distribution to be y.i/2{Ti — 0.5)Af G (-54,-20) ms. 
Thus, the WMAP 5t = 0.5 hypothesis can be conser- 
vatively rejected at 99.999% confidence. The same non- 
parametric method gives the (2.5%, 97.5%) confidence in- 
terval to be (—47, —26) ms, i.e. taking into account Bennett 
et al. (j2003ap 's (Sect. 6.1.2) estimate of the uncertainty in 
the timing as 1.7 ms, the hypothesis of a —25.6 ms offset is 
consistent with the data. 



3.1. Quadrupole overestimation induced by the calibration 
error 

Unfortunately, the simple procedure presented here is only 
sufficient to provide evidence for the timing-offset-induced 
calibration error, not for compiling improved maps directly 
from the calibrated TOD. As explained in Sect. 12.11 a 
timing-offset-induced calibration error on a signal domi- 
nated by the dipole can be approximately reversed, since 
the calibration is itself based on a model dipole that is a 
good approximation to the measured dipole. However, the 
remaining signal is not a dipole, so it is more difficult to 
correct the remaining components of the wrongly-calibrated 
signal. As shown in Roukema (j2010|) . the present method 
of reversing the calibration error, i.e. the use of g~^, in- 
troduces a blurring effect on small (^ 4') scales. To obtain 
CMB sky maps with the calibration error fully removed, it 
is not obvious that there is a simpler alternative to redoing 
the full calibration of the uncalibratcd TOD. 

Nevertheless, an approximate reversal of the —25.6 ms 
timing-offset-induced calibration error in at least the dipole 
component of the signal is possible using the same data 
analysis pipeline as above. In this case, by how much does 
the quadrupole decrease when the artificial quadrupole-like 
signal is approximately removed? 

Using the presently available version of the WMAP 
3-year, calibrated, filtered TOD in the W band, i.e. the 



band where the cosmological signal is the least affected 
by foregrounds, the quadrupole can be calculated consis- 
tently for maps at both 5t ~ 0.5, the incorrect timing 
offset, and 6t — 0, the timing offset that yields approx- 
imately corrected maps. Using maps made at a resolu- 
tion of A^sidc = 512, and subtracting the 3-year maxi- 
mum entropy synchrotron, free-free and dust maps pro- 
vided by the WMAP coUaboratior FH. p seudo-C/ estimates 
(after monopole and dipole removalri from the 3 years, 
for the 4 W band DA's, are (3/7r)C2 = 170 ±AfiK^ and 
103 ± 4^K^, for 6t ~ 0.5 and St = 0, respectively, for the 
KQ85 sky mask, and {3/tt)C2 = 79 ± 3 //K^ and 41 ± 2 yuK^ 
for the KQ75 sky mask, where errors are standard errors in 
the mean over the 12 DA/year combinations, and 80 map 
iterations are made for each DA/year combination. The 
dipole removal has little effect on these estimates. Removal 
of only the monopole gives (3/7r)C2 = 170 ± 4/iK^ and 
103 ± 4^K^, for St ~ 0.5 and St = 0, respectively, for the 
KQ85 sky mask, and (3/7r)C2 = 78 ± 3 //K^ and 40 ± 2 /xK^ 
for the KQ75 sky mask. Hence, the —25.6 ms timing-offset- 
induccd calibration error implies that estimates of the CMB 
quadrupole amplitude (3/7r)C2 based on WMAP 3-year W 
band sky maps derived from the incorrectly calibrated TOD 
are overestimated by 64 ± 6% and 94 ± 10% for the KQ85 
and KQ75 sky masks, respectively (for removal of both the 
monopole and dipole). 

This is a significant and substantial systematic error. 
Given the differences in method and data subsets [e.g., Liu 
& Li (|2010ap estimate cross-quadrupoles of the V and W 
five-year data, while the present estimate is the W 3-year 
auto-quadrupolc] , the KQ75 drop in quadrupole power ap- 
pears to be roughly consistent with that of Sect. 5.3 of Liu 
& Li (|2010ap . i.e. greater than that of Moss et al. (P(JTU| 's 
toy model estimate (Fig. 1 caption). Moreover, the problem 
of large-scale CMB power being concentrated towards the 
Galactic Plane is exacerbated with this approximate cor- 
rection of the calibration error. The Galactic Plane is the 
part of the sky where difficult-to-estimate systematic error 
can most be suspected. 



4. Discussion 

4.1. Sensitivity to fitting method 

These results do not appear to be sensitive to the choice of 
symmetrical fitting function. Least-squares hyperbolic best 
fits, i.e. best-fits oi ai{St)^ + biSt + Ci to [si{St) + 0.5\'^, where 
the vertical offset of 0.5 and positive square root of the 
best fit give a top quadrant hyperbola, yield similar values 
to those in Table [D For example, over all wavebands, this 
hyperbolic fit gives ^1/2 (t;) = 0.22 ± 0.04 and /ii/2(Ti — 
0.5)At = -29.4 ± 4.3 ms. 

Since the St values used in sampling are symmetric 
around St = 0.5, it is difficult to see how this could induce a 
bias against the St = 0.5 hypothesis. Could it lead to an in- 
accurate estimate in the case that (r^ — 0.5)At = —25.6 ms 
is correct? Table [3] shows that using a St sample symmet- 
ric around St = gives results compatible with those in 
Tabled 



^^ http://lainbda .gsf c. iiasa.gov/prodvict/map/dr2/ 1 
mem_maps _get . cf m 

^^ The ANAFAST routine of the GPL version of the HEALPlx 
software (Gorski et al. I1999P was used here. 
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Table 3. Quadratic fit estimates [median ^i/2(Ti — 0.5)Ai and 
and standard error in the median, in ms] of minimum standard 
deviation, as for Table [J] using 5t intervals symmetric around 
St = 0.0 in the range —5 < 5t < 5, i.e. possibly showing bias 
against St = 0.5. 



units : 
band 


^1/2(1-0 0-Ml/2(^>) 

At At 


At fii/2(n ' 
ms 


-0.5)At cr^^^^(,^)At 
ms ms 



K, Ka 

Q 

V 



0.20 

0.21 

0.15 

-0.17 



0.07 128.0 

0.04 104.2 

0.13 76.8 

11.56 52.1 



-37.9 

-30.3 
-26.7 
-35.0 



9.2 

3.7 

9.8 

602.3 



all 



0.21 



0.03 



-30.8 



3.2 



^In W, the best-fit parabola for one of the year/DA combina- 
tions had a maximum and was ignored. 
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Fig. 5. Mean and uncertainty of Si{5t) in the W band including 
a posteriori extra maps for large values of \5t\, as for Fig. [2] 



Could the dependence of Si on 6t be asymmetric? The 
toy model (Sect. [121) does not suggest any obvious asymme- 
try, and it is not obvious vifhat detailed differences between 
the toy model and the real data could be sufficiently asym- 
metric to perturb the method. Moreover, the WMAP scan 
paths are curved (non-geodesic) paths on the 2-sphere that 
have many symmetries. Nevertheless, it is not obvious that 
Si should vary with an exact symmetry around its mini- 
mum, even if approximate symmetry seems reasonable. If 
the dependence of Sj on St approached linearity for large 
\St\, then one way of testing this could be to make linear 
fits to subsets of data for large \St\. However, comparison of 
parabolic and hyperbolic fits favours the former, i.e. Si(St) 
does not appear to approach linearity for large |i5i|. Thus, 
linear fits for large \St\ would depend strongly on the 6t do- 
mains chosen. If an asymmetric model for Si dependence on 
St could be found that is close enough to parabolic in shape 
in order to provide fits in K and Ka, Q, and V as good as 
those shown in Fig.[2j then it could conceivably be possible 
that the St = 0.5 hypothesis could be made consistent with 
the resulting minima. However, in that case, the fact that 
the symmetry assumption gives a solution consistent with 
the timing offset recorded in the TOD files would have to 
be attributed to coincidence. 



Table 4. Quadratic fit estimates for usage of the calibrated 
TOD with an extra term, i.e. using the 80-th map iteration, as 
for Table [21 for the bands with shorter calculation times. 
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0.06 
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Q 




0.40 


0.13 


104.2 






-10.0 


13.7 


V 




0.11 


0.65 


76.8 






-30.3 


49.7 


all 




0.33 


0.05 


- 






-20.5 


6.4 



4.2. Consistency between wavebands 

In Fig. [21 the median standard deviation s, is most noisy in 
the W band and a little noisy in the V band. The quanti- 
tative best fit minima shown in Table [D significantly favour 
the —25.6 ms offset hypothesis and reject the other hy- 
potheses. The choice of intervals for mapmaking was opti- 
mised for testing the St = and St = 0.5 hypotheses rather 
than an offset in ms, which is why the spreads of StST for 
V and W are less than for K, Ka and Q. Thus, with hind- 
sight, given the results found here, a possible explanation 
for the higher noisiness in V and W would be that calcula- 
tions were not spread over a wide enough interval in St in 
these bands. 

Another consistency check between results in different 
wavebands is that if the increase in Si away from a preferred 
timing offset is mainly an effect of incorrect dipole-bascd 
calibration of the observed dipole signal, then this depen- 
dence should only be weakly dependent on waveband. 

Figure [D shows the same information as in Fig. [U with 
zeropoint removal but no scaling of the standard deviations, 
i.e. CTi — Ai, against timing offsets in milliseconds centred 
on the hypothesis preferred by the data, i.e. {St — 0.5) Ai = 
—25.6 ms. In the range of approximately {St — 0.5)At £ 
(—25.6 ± 175) ms, the V and W band results appear fully 
consistent with the K and Ka, and Q results. The V and 
W bands do not appear to be more noisy than the K and 
Ka, and Q bands. Moreover, the dependence of the effect 
on {St — 0.5)Ai appears to be approximately independent 
of waveband, at least in this central range. A posteriori 
calculation of maps in W for larger absolute values of St 
(e.g. to about ±640 ms as in the K and Ka bands) could 
reasonably be expected to lead to a smaller uncertainty in 
this band. 

This is indeed the case. Extra maps were calculated for 
St ~ — 10, . . . , —6, 6, . . . , 10 in the twelve W band DA/year 
combinations. If all the W maps are analysed together, 
the resulting best estimates are iii/2{Ti) = 0.19 ± 2.24 and 
Mi/2('''i — 0.5)At = — 16.1±116.8 ms. However, Fig.[5|shows 
that the dependence of si on St in the W band is more 
complex than the quadratic dependence that gives visually 
acceptable fits in the K and Ka, Q, and V bands. 

Figure[4|also shows that beyond the central range where 
the different bands appear consistent, there appear to be 
significant differences in the amplitude of the variation in 
Si between the different wavebands, given the uncertainties 
estimated among DA/year combinations within each wave- 
band, even though the estimates of minima are clearly con- 
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sistent (Table [5]). This variation among wavebands might 
be due to an increased role of non-dipole components of the 
signal, e.g. Galactic foregrounds that are brightest in the 
lower frequencies. Since cosmological perturbations should 
be much weaker than both the dipole and low frequency 
foregrounds, the approximate frequency independence of 
the effect near the minimum variance and frequency de- 
pendence away from the minimum is qualitatively consis- 
tent with the effect near the minimum being dominated by 
the dipole, and the effect further from the minimum in- 
volving coupling effects between the dipole calibration, the 
measured dipole, and foregrounds. The beam size and side- 
lobe effects reported by Sawangwit & Shanks (j2010[) and 
Liu & Li (|2010b[) may be related to this latter speculation. 

4.3. Sidelobe and other beam effects 



Sawangwit & Shanks (|2010p have estimated that WMAP 
beam profiles at faint flux density levels, especially in the 
W band and weakly in the V band, are wider than those 
estimated by the WMAP collaboration for bright flux den- 
sity levels. Liu & Li (|2010bp have suggested that WMAP 
sidelobe interactions with the dipole signal create an effect 
that may be effectively modelled as if it were a positional 
offset. From Fig. 1 of Page et al. (|2003)) and Fig. 2 of Hill 
et al. (Unnni), for the first year and 5-year WMAP data re- 
leases respectively, it is clear that sidelobes in the beam 
shapes are strong for the W DA's, weak for the V DA's, 
and much weaker for the K, Ka, and Q DA's. 

The beam shape dependence on waveband suggests an 
explanation not only for the variation in si dependence 
on {5t — 0.5)Ai among wavebands at large \{5t — 0.5)Ai|, 
but also for the complexity in the shape of Si{6t) in the 
W band (Fig. [5]) and the associated large uncertainty in 
A«i/2(t-j - 0.5)At in the W band (Table ^. The (partial) 
invertibility of a timing-offset-induced calibration as pre- 
sented in this work depends on the details of the scanning 
pattern. For a complicated beam profile, the relation be- 
tween the sidelobe orientation, the scanning pattern, and a 
slightly displaced dipole, is likely to be complicated. For a 
simple profile, the relation is likely to be simpler, i.e. less 
noisy. Thus, the complexity in Si{dt) in the W band is qual- 
itatively consistent with Liu & Li (j2010bP 's estimate that 
sidelobe interaction with the dipole is important in map- 
making, given the WMAP collaboration's estimates of the 
beam profiles shown in Fig. 1 of Page et al. (j2003[) and 
Fig. 2 of Hill et al. (|2009p . It is also qualitatively consis- 
tent with the result of Sawangwit & Shanks (|2010p that 
the difference between their radio-source-based beam pro- 
file estimates and the WMAP collaboration's Jupiter-based 
estimates is strongest in the W band. 

4.4. Extra term from iterated maps 

Although use of an iterated map is equivalent to adding an 
extra term to Eqs (J13p and (jl4[) . introducing more noise into 
the method, the result should be statistically compatible 
with the above results. Checking the bands with fewer num- 
bers of DA's and observations, i.e. adding the 80-th map 
iteration to these equations, took about a calendar month 
of calculations. The results are shown in Table HI These 
mainly consist of an increase in the noise, with a slight in- 
crease in 6t values, i.e. ^i/2(Ti — 0.5)Ai = —20 ±6 ms. This 



estimate still disagrees with 6t = to high significance, and 
agrees with (6t — Q.5)At = —25.6 ms at O.Scr. 

4.5. Devil's advocate: could the calibration be correct? 

It is difficult to see how the variance per pixel can be de- 
creased by an arbitrary shuffling of the TOD. The simplest 
explanation for the best estimate of yUi/2(''"i — 0.5) At = 
—29.7 ± 2.5 ms is that the —25.6 ms timing offset correc- 
tion was not implemented when the WMAP uncalibrated 
TOD were calibrated. However, could it be possible to avoid 
this inference, i.e. could the calibration be correct? Some 
possibilities and counterarguments arc as follows. 

(i) As mentioned above, a model in which the dependence 
of Si on St is asymmetric might, in principle, give the 
minimum variance per pixel to occur near Ti = 0.5. This 
would require a heuristic model for the asymmetric be- 
haviour, and at least qualitatively would seem to require 
a much more complex shape than a parabola. 
(i.l) In this case, it would be a coincidence that the 
symmetry assumption gives a result consistent with 
one of the two versions of Liu et al. (|2010ap 's hy- 
pothesis, as listed in Table [TJ That is, the successful 
prediction of one of the two versions of Liu et al. 
(|2010a|) 's hypothesis would be a coincidence, 
(ii) A heuristic model could relate to the fact that the K, 
Ka, and Q bands have much stronger Galactic con- 
tamination than the V (and W) bands. Could the 
l^i/2{Ti — 0.5) At = —29.7 ±2.5 ms estimate be primarily 
an effect of Galactic contamination? 
(ii.l) The dependence of Galactic contamination on 
waveband is very strong going from K to W. Tablcs[2] 
and Tabic [3] do not show any obvious trend in this 
sense. 
(ii.2) The statistic Si is a normalisation of Ci, which is 
the median^ over the pixels of a given map, of the 
standard deviation per pixel. For Galactic contam- 
ination to affect the median, by shifting flux den- 
sity onto or off given pixels, it would need to affect 
about half of the pixels at high Galactic latitude. 
However, a —25.6 ms timing offset corresponds to 
only 4', while pixel sizes are about 7.3°, i.e. about 
four orders of magnitude greater in solid angle. Even 
in the K and Ka bands, where the Galactic contam- 
ination is very strong, the number of tiny positional 
shifts of flux density onto or off a given pixel should 
be very high, so that the median over all valid pixels 
could be expected to vary smoothly and symmetri- 
cally around the correct timing offset. 
(ii.3) The coincidence argument also applies here. The 
Galactic contamination would have to mimic one of 
the two versions of Liu et al. (|2010ap 's hypothesis 
by chance, in sign (t; < 0), in amplitude, and in 
having an approximately constant value expressed 

as ^il/2{Ti - 0.5)At. 

(iii) Could the small-scale blurring effect (Roukema 120101 
and references therein) have played a significant role in 
this analysis? 

(iii.l) If it did, then it would have caused the results to 
be biased in favour of r^ = 0.5, i.e. ^i/2('''i — 0.5)At = 
ms. In other words, correcting the present result 
for the small-scale blurring effect would shift it in the 
opposite direction, giving /ii/2(Ti — 0.5)At < — 29.7± 
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2.5 ms. This would worsen the rejection of the tj = 
0.5 hypothesis, 
(iv) Liu & Li (|2010ap showed and Moss et al. (PUTU)) con- 
firmed that the CMB quadrupole is reduced in ampli- 
tude for At{6t — 0.5) ~ —25.6 ms. If this is just a coin- 
cidence that relates the CMB quadrupole, the velocity 
dipole and the WMAP spacecraft scan pattern, could 
the variance minimisation method be strongly biased 
towards minimisation of the quadrupole, so that these 
constitute two instances of a single coincidence? For ex- 
ample, let us suppose that Ji = 0.5 is the correct value. 
It is already known that in this case, the mean signal 
at 5t = will be smaller. Could this imply that the 
standard deviation per pixel at 6t = would be smaller 
too? 

(iv.l) If dt = 0.5 were really correct, then, although 
the mean signal at 5t ~ will be smaller, this 
would only be a result of the fortuitous cancelling, 
in the mean, of a component of the real signal by the 
dipole-induced difference signal. However, reducing 
the mean by the addition of anticorrelated errors 
does not reduce the standard deviation per pixel. 
The scatter induced by a wrong assumed value of St 
is not cancelled by a reduced mean value. 
(iv.2) Calibration errors induced by a wrong value of 
St apply to the full CMB signal, not just to the 
quadrupole part of the signal. The latter is a very 
weak part of the full CMB signal. Thus, scatter in- 
duced by a wrong value of St can have only a weak 
dependence on the quadrupole. Table [2] shows that 
St = 0.5 is still strongly rejected in the K and 
Ka, and Q bands, where galactic foregrounds play 
a strong role. 

Retaining the hypothesis that the calibration was cor- 
rect does not seem to be easy. 



5. Conclusion 

While the —25.6 ms offset between the times in the Mcta 
Data Set and the Science Dat a Set in the WMAP TOD 
files, discovered by Liu et al. (|2010ap . did not lead to an 
error in compiling the calibrated, filtered, 3-year TOD into 
maps (Roukema l2010p . it is difficult to avoid the conclusion 
that this timing offset did induce an error at the calibra- 
tion step. The two different versions of Liu et al. (|2010ap 's 
hypothesis gave numerical predictions, listed in Table [1] 
The constant offset in ms hypothesis is consistent with the 
analysis of the variance maps, summarised in Table [2] This 
result does not seem to be sensitive to the fitting method 
(Sect. HTTj) . it is consistent among all wavebands (Sect. H?^ . 
using highly iterated maps adds noise and gives consistent 
results (Sect. 14. 4p . and alternative explanations seem to 
be speculative and/or fail to save the St = 0.5 hypothesis 
(Sect. HTS]) . The timing offset that is numerically implicit in 
the TOD files (e.g., see Appendix A, Roukema l2010[) would 
appear to provide the simplest explanation. A new analysis 
of 7 yea rs of Q, V, and W band WMAP TOD by Liu ct al. 
(j2010bp finds a similar result. Thus, it appears that (1) in 
the calibration step, both the spacecraft attitude quater- 
nion timestamps and the observational timestamps were 
used, for the dipole and observational data, respectively, in- 
ducing a calibration error, while (2) in the mapmaking step, 
the observational timestamps were (correctly) assumed to 



be correct both for dipole and observational data, inducing 
no further error, but retaining the original calibration error. 
This supports the claims by Liu et al. (Liu & Li l2010al Liu 
et al. I2010ap that the CMB quadrupole has been substan- 
tially overestimated. A rough estimate made here is that the 
quadrupole estimated in maps based on the wrongly cali- 
brated TOD is overestimated by about 64 ±6% to 94 ±10% 
for the KQ85 and KQ75 sky masks respectively (Sect. [XT]) . 
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